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1. Introduction 

Random gaussian fields constitute an important paradigm in physics. They serve as 
a simplified model for diverse physical ingredients such as the energy landscape of 
glassy systems [1] [51 [31 H] , wave functions of quantum systems with chaotic classical 
dynamics [S], and the fluctuations of the Cosmic Microwave Background radiation 
[6]. The critical points of these gaussian fields - minima, maxima and saddle points 
- characterize the behavior of the entire field and thus carry a large amount of 
information. For instance, large scale structure in the universe can be analyzed using 
the critical-point density [7]. In other cases critical points are directly related to 
physical phenomena such as reflection patterns of sunlight from the sea surface [S], 
and the onset of glassy behavior in a complex energy landscape, see [1] and references 
therein. An additional novel context where critical points play an important role is 
that of cold atoms in disordered optical potentials [9l[T0l[Tl]. In these systems, proper 
detuning of the light frequency forces the atoms to gravitate to the local extremum 
points of the intensity, so that statistics of the critical points map directly to statistics 
of the atom positions. 

Most studies of critical points have focused on their density. However, in order 
to characterize their statistical properties it is also important to understand their 
correlation functions. In particular, in order to calculate fluctuations of the number of 
critical points, within a given volume of the system, one needs the two-point correlation 
function. 

In this work we shall address the problem of pair correlations of critical points. 
The one-dimensional version of this problem was solved many years ago |12| I13| . In 
two dimensions, some asymptotic expressions are known |14| I15| . along with some 
numerical data [16] (see also jl7| I18|). For a comprehensive review, concentrating on 
wave systems, see [12] and references therein. We shall provide a detailed derivation 
for the general asymptotic formula of critical point correlations in any dimension 
d > 2, for both long and short distances. Some of this work, especially in the long- 
distance limit, may be seen as an expansion and generalization of methods developed 
previously, see e.g. [20 l [2T 1 [T4] . 

There are several types of densities that can be associated with critical points. 
In order to present them, let 0(r) be a gaussian random fleld in d spatial dimensions, 
r = (ri, r2 • ■ ■ r^), with zero mean and an isotopic correlation function 

(0(r))=O, (<^(r)<^(r')>=G(|r-r'|), (1) 

We shall assume that the second derivative of G{r) is finite at r = and vanishes at 
r — >■ oo. We note that in general this implies that r'^^G^r) = vanishes in the limit 
r — >■ cx), for < 7 < 2. 

Critical points are the points where the gradient of the fleld vanishes, and their 
density is given by the well known Kac-Rice formula: 

p„,(r) = |detiJ|,5(V</)) (2) 

where H is the Hessian matrix whose elements are the second derivatives of the fleld: 

This density, which we shall refer to as the unsigned density of critical points gives 
an equal weight and sign to each one of the critical points. There are, however, many 
other possible choices of density. For instance one may define the signed density as 
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in Eq. ^ but omitting the absolute value operation. In two spatial dimensions this 
density assigns a plus sign to minima and maxima points and a minus sign to saddle 
points. Another possibility is to define a density of only minima points of the gaussian 
field. The main focus of this work is the correlation functions of the unsigned density 
and the density of minima points: 

C(r) = (p(r)p(O)) (4) 

where p(r) denotes one of these critical-points densities. However, the methodology 
we shall develop in order to calculate these functions can be generalized in a 
straightforward manner to the correlation functions of different types of critical point 
densities. 

The central result of this work is that the asymptotic behavior of the correlation 
function is 

^^"^'-X a3r2-^+'= 0<r«a 

where {p) is the average density of critical points, ai,a2 and are constants which 
depend on the type of the critical-point density and the dimensionality of the system 
d, 

d^Gjr) 
dridrj 

is the Hessian matrix of the correlation function of the gaussian field, and a is its typical 
length scale. In this formula fc = corresponds to the correlation of the unsigned 
density, while fc = 3 corresponds to the correlation of mimina points. Notice that the 
correlation function also contains a (5- function at the origin, C (r — > 0) = {p)d (r), due 
to trivial self correlations of the critical points. 

The paper is arranged as follows: In Sec. [2] we shall sketch the derivation 
of asymptotic form of the correlation function ^ deferring details to [Appendix A[ 
Then we shall apply these results to find the correlation functions in two and three 
dimensions, in the long-range asymptotic limit. Finally, we will discuss some examples 
and conclusions in Sec. |4l Appendix A contains the technical details of the derivation, 
and Appendix B provides information regarding the numerical calculations which were 
performed in order to compare Eq. ^ with the exact result. 



lHG(r)].,, = (6) 



2. Asymptotic expressions for C{r) 

In this section we shall sketch the calculation of the asymptotic behavior of the pair 
correlation function of critical points. We start with a general formalism in which 
we define the various ingredients of the problem needed for our calculation. Then we 
show how to obtain the long range and the short range asymptotic behavior of the 
correlation function of critical points. 



2.1. General formalism 

Consider the critical points of a d— dimensional homogenous and isotropic Gaussian 
random field (/i(r), r G E"*. The unsigned critical point density is defined by Eqs. 
([3]) and Here the Dirac 5- function (5 ( V0) selects points with zero gradient, and 
the Hessian determinant provides the correct normalization since the integral over an 
infinitesimal region around the critical point is unity. 
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This density does not distinguish between extrema and saddle points. However, 
in this work we shall also be interested in the density of the local minima of the field 
only. This density may be expressed as 

Prmn (r) = det i7 [] 9 (dct 5 iVcj,) , (7) 

3=0 

where dct _ff is the j-th minor of H, i.e. the determinant of the matrix obtained 
when the first j rows and columns of H are removed (the zeroth minor is defined as 
= /f). To understand this result notice that minimum points are associated with 
the case where all eigenvalues of the Hessian matrix Xj [j = 1, • • • , d)) are positive 
(Note that the Hessian matrix is symmetric and therefore its eigenvalues are real). 
But this condition can be expressed in terms of determinants of the minors since 

n -=1 e (A,) = n,ti e (n a,) = n,to © {detHU)). 

The translational invariance and isotropy of the correlation function of the 
gaussian field implies that the pair correlation function of critical points depends only 
on the distance between the two points and therefore one may arbitrarily choose these 
points to be: one at the origin, ri — (0, • • • , 0), and the other at r2 = (r, 0, • • • , 0). 
Then in order to perform the average the product of two densities over the gaussian 
field it is sufficient to consider the joint probability distribution of the gradients of the 
fields and the second derivatives at the two previously chosen points ri and r2. It is 
thus convenient to define the vector: 

* = (V(/)(ri),tAi,V0(r2),iA2)): (8) 
where xp^ is a vector constructed from the elements of the Hessian matrix, i.e. the 
second derivatives of the gaussian field: 

tA,. = ({iJ.,,(r.)}) (z<j). (9) 

Being expressed in terms of derivatives of a Gaussian field, the vector ^ 
is normally distributed, and thus defined by its mean (which vanishes) and the 
correlation matrix: 

B=(*®*) (10) 

where (g) denotes an external product. Thus 

where 2n ~ 2d + {d — l)d is the number of components of the vector and d"^ 
represents the an infinitesimal volume in the 2n-dimensional space. 
Now the correlation matrix has the general form 

f 6(0) b{r) \ 

where the matrix elements of 6(0) (which is independent of r) correlate components 
of 4' at a single point in space, while those of 6(r) (and its transposed matrix 6*(r)) 
correlate components of the field associated points that are a distance r apart. This 
structure is convenient for an asymptotic expansion since at distance greater than the 
correlation length of the gaussian field 6(r) — > , while in the opposite limit, r — >■ 0, 
the matrix 6(r) converges to 6(0). Thus 6(r) and 6(r) — 6(0), respectively, may be 
treated perturbatively in order to obtain the asymptotic behavior of C{r). 
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In order to take advantage of this property it is convenient to use the Hubbard- 
Stratonovich transformation and represent the integral (|lip in the form: 

C{r) = I d*p(0)p(r) J dUcxp (-^nSH + ■ .(13) 

The integration over the components of ^ associated with the field gradients 
can be readily carried out, and by performing also the integration over their Fourier 
conjugates we obtain: 

^^'^ = (2^)^»-VdetA / / ^"'^"P (-r^" + • ^) (^'^ 

where p(r) — p(r)/5(V0(r)), and the remaining integral is over the vector 

V' = (Vi,V'2) = (V'(ri),V(r2)), (15) 

and its Fourier conjugates which we denote by tt. The matrices A and B may be 
expressed in terms of submatrices of B. Details on their structure may be found in 
[Appendix A[ 

2.2. Long range asymptotics 

In order to obtain the asymptotics in this limit we use the fact that the off-diagonal 
blocks of B approach in the limit r — > oo, and therefore may be treated as small 
parameters. It is thus convenient to introduce a dummy parameter to B |20| : 



Tjb\r) 6(0) 

T] will serve as a bookkeeping parameter of the pertubation theory, and at the end of 
the calculation will be set to 1. With this addition, A and B appear in the form 

A-( \ / a{r)\ 

^ ~ 1^ a(0) )^'\ a{r) J ' 

[ 6(0) r 'UM j H d2ir) d,{r) j ^^"^ 

where a(r), 6(r), di{r), and d2{r) are submatrices which approach zero as r — > cxj (see 
Appendix A for details). Substituting these matrices in and expanding in rj give 
the formal expansion of the correlation function, 

C{r) = Co + 7,Ci(r) + ifC^ir) + • • • , (19) 

where 

(p(G)p(r))o 2 
^" = (27r)'^deta(0) = ^'^^^ 
is the disconnected part of the correlation function, and here we use the notation 

1 
2 



0= , ^ _ = I dV(---)expf-^i/>^o'V) , (21) 



dei(27rBo) 
where Bq = -B|,,=o- 

The leading order corrections are obtained by expanding the integrand of (jl4p in 
?7, and performing the resulting integral. The result takes the form 

Ci{r)^\i:r{KF), (22) 
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where 

b{r) 



^ = 4-M ,L T ^-^ (23) 



and 



The leading contribution to C2 depends on the behaviour of G{r) and its derivatives. 
If, for simphcity's sake, we assume that G{r) does not have a characteristic wavelength, 
i.e is not oscillatory or exponentially decaying, then we have: 

C2{r)^ColTr[aiO)-'a{r)]\ (25) 

If G(r) has a characteristic wavelength a more complicated expression is necessary, 
and appears in (jA.lOp of the Appendix. 

This result gives the leading order behavior of C(r) in the long range asymptotic 
limit. Further details on the structure of K, F, and the matrices a(r) and b(r), are 
left for [Appendix A| For now we just point out that a{r) = —Hair), while V^G(r) 
can be constructed from elements of b{r). Thus performing the trace operations in 
Eqs. ([211), (HH) leads to our main result, Eq. ([Sj), with the constants: 

«i=Co/^[V^G(r)|^J-^ (26) 

a2 - Cod [TvHG{0)^r' 

where / vanishes for the correlation function of the unsigned density, while it is a 
constant of order unity for the minima correlation function. This constant, which 
depends on the the dimensionality of the system, will be discussed in Sec. [3l below. 



2.3. Short range asymptotics 

As r — > 0, the matrix B becomes singular. As a result, both A and B in Eq. ([T^ 
become singular as well. Therefore, in order to calculate the correlation function, one 
should handle carefully the divergency of the prefactor in Eq. (fT4|) in the limit r — s- 0, 
and cast _B in a form appropriate for expansion. This process is rather technical, and 
we delegate the body of this discussion to [Appendix A[ bringing a bare outline here. 
We choose r itself to be our expansion parameter. It can then be shown that for 

[dct A]-^^'^ (X r-"^ + ■ ■ ■ (27) 

Now, as r — >■ 0, the matrix difference b{r) — b{0) approaches zero. Therefore, in this 
limit it is convenient to represent the matrix B in the form: 



B = 



61 (r) &i(r) 
6i(r) fei(r) 



Ab 
Ab 



with bi{r) - 
implies Ab 



6(0)-di(r),A6 
> when r 



b{r) - 6(0) - [di{r) - d2{r)] (see Eq. 
0. We then perform the unitary rotation 



(28) 



13). This 



V2 V 1 -1 



i>2 



(29) 
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and a similar rotation to the Fourier conjugate variables tt. This transformation block- 
diagonalizes the matrix B, so that C{r) takes the form but with the matrix B 
replaced by 

An additional transformation block-diagonalizes bi + Ab, and allows one to cast the 
resulting matrix as a power series in r and perform the integral. The result is 
w f ''^ unsigned density 

" I r° density oi mmmia ^ ' 

Combining Eqs. (|14|) . ([27| and ((3T|) leads immediately to our main result, Eq. ([5]). 



3. Application in two and three dimensions 

In this section we shall find an asymptotic expression for the correlation in two and 
three dimensions. We begin by calculating Cq and F given by Eqs. (|20p and 
respectively. Both formulae no longer possess an explicit dependence on r. As a 
result, the isotropy of the field greatly simplifies calculations. Thus, finding Co 
clearly reduces to simply finding the density expectation: 

Co (X (p(0)p(r))o = {m)l (32) 
Isotropy also implies that all non-zero elements of F appearing in Ci have the 
same value /i, thus: 

/i = co/cx ((m^'' 



drf 



r=0/ 



= l((p(0)Trif(0))o)' (33) 

where / is the constant which has been introduced in Eq. (j26p . and Fi^n-d+i is an 
example of one of the relevant non-zero elements of F. 

The density (p) (and thus Co) has been calculated by others [HIT], therefore, it 
remains to find /. We devote the rest of this section to this calculation, and to be 
concrete we choose the correlation function of the gaussian field to be 

d-2 

Gdir) - (^) ' Ji^ (vrr) , (34) 



2r^ 

where c? = 2, 3 is the dimensionality and is the Bessel function of order v. This 
form of correlation appears naturally in the context of ballistic chaotic wave systems, 
the so-called Random Wave Model [SI [IH] ■ This choice does not affect the calculation 
of / but is useful for comparison with numerical results. 



3.1. Two dimensions, d = 2 

In two dimensions ij) = {(jy^x, 'Pxy, fpyy)-: where for brevity we use x, y subscripts to 
denote derivatives with respect to the corresponding coordinates. In this case we 
have defined xjj = xj^i since ■02 no longer appears in expressions and ([55)). By 
transforming to the normalized variables [8l I22|. 

X (X (j)xy 

Z CX (jixx + (Pyy 
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one obtains dcti? oc 2Z^ — — Y^, Tr H (x Z, and the Gaussian measure becomes 
diagonalized. It is now possible to calculate {p{0) Ti H (0)) ^ by transforming to 
spherical coordinates. A similar derivation is found in |15j and we omit the details. 
The final result is / = for the correlation function of unsigned density. (This is in 
fact true for any dimension and can be seen from symmetry considerations.) For the 
minima-minima correlation, 

We conclude this section by referring to figure [T] This shows a comparison of our 
analytical asymptotic expression, and a numerical evaluation. Details of the numerical 
procedure used in order to get this figure can be found in Appendix B. 

3.2. Three dimensions, d = 3 

The critical-point density was found for d = 3 in the context of the Cosmic Microwave 
Background radiation [7]. In this case there are six variables in the Hessian: 
tp = {4'xx, (f'xy, (^xz, 't'yz, (f'zz), where as in the two-dimensional case, we write for 
brevity xp = xpi and the subscripts denote partial derivatives. Here we transform the 
Gaussian integral over ^p into an integral over the eigenvalues of the Hessian matrix 
H. To this end it is convenient to represent Gaussian exponent in the form [HIT]: 



^^'^ 4(02^) V d + 2 

where here we have used the block form of Bq in Eq. psp . Employing the standard 
approach of random matrix theory |23| one may transform from the integration 
variables ip to to an integration over eigenvalues oi H, \ = (Ai, A2, A3). The measure 
of the integral in this case is 

dVexp(^V[6(0)]"^-i/') = (36) 
n ■=! d\, 



(A2 - Ai) (A3 - A2) (A3 - Ai) exp 



with Mij = Sij — i. Following Ref. [7], further diagonalization of the Gaussian 
exponent followed by a change to spherical coordinates allows one to calculate 
(p(0) Tt H{0)) and to deduce the value of the constant / from Eq. ([55)1 . The result is 
/ = for the correlation function of unsigned density, while for the minima-minima 
correlation: 

/ - ''''' (37) 
9(29-6^6)277 

We conclude this section by referring to figure [5J which shows a comparison of 
our analytical asymptotic expression, and a numerical evaluation. 

4. Discussion 

In this work we calculated the asymptotic behavior of the pair correlation function 
of critical points, C(r) in two limiting regions, r — > 00, and r — > 0. We have 
shown that in both limits the behavior correlation function depends on the types of 
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(a) Minima-minima correlation function 
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Figure 1: The critical point correlation functions for a gaussian random field in two 
dimensions, with Bessel function correlation G (r) = Jq (irr). The term Cq represents 
the disconnect part of the correlation function, while Ci and C2 are the leading terms 
of the long range asymptotic expansion. The black dots represent the exact result 
obtained by numerical calculation explained in Appendix B. 



densities of critical points, i.e. whether it is the unsigned density or the the minima- 
minima density. In the short range asymptotic it is manifested through the power law 
dependence of the correlation function: C{r) cx r'^~'^ for of the unsigned correlation 
and C(r) cx r^^'^ for minima-minima correlations. 

In the long range asymptotics we have seen that there are two terms contributing 
to the correlation functions, see Eq. ([5]). However, for the unsigned correlation the 
prefactor of the term proportional to V^G'(r) vanishes while for the minima-minima 
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(b) Unsigned correlation function 

Figure 2: The critical point correlation functions for a gaussian random field in three 
dimensions, with G(r) = SiSJii^ xhe meaning of lines and the dots in this figure is 
the same as in Fig. 1 



correlation both terms contribute. 

Let us consider some specific examples in order to clarify the nature of the 
contributions of the two terms, V^G{r) and Tr [_ffG(r)]^, to the long range asymptotic 
limit. Consider the case where the correlation function of the gaussian field is an 
oscillatory function of the form 

hm Gir) A^-^ (38) 
1 — >oo r' 

where A is the amplitude, cr is a constant which sets the length scale of the oscillations, 

while 7 > determines the behavior of the envelope of the oscillations. This type of 
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correlation describes, e.g., a field generated by wave chaos For this case in the 
asymptotic limit ?' — > oo 

Acj^ (39) 

while 

Tr[i7c(.r-^^^ (40) 

Thus for large enough r the first term, Eq. (|39p . dominate the behavior (unless its 
prefactor ai vanishes as for the unsigned correlator). 

Consider the case where at large distances the correlation function is described 
by a simple power law behavior, i.e. 

lim G{r) ^ (41) 

r— >-oo r> 

where 7 > —2. This type of correlation, in particular for 7 < 0, is believed to describe 
the potential landscape of glassy systems [53]. For this case we find that 

V4G(r) ^ 7(2 + l){d - 7 - m - 7 - 4)-^, (42) 

while 

Thus the term is dominant for 7 > 0, while is dominant for 7 < 0. 
Finally, for logarithmic correlation: 

lim G{r) — > Alog(r) (44) 



Tr[i7G(r)]'^7'(7 + rf)'3^- (43) 



we have, 



while 



V4G(r) ~ 2(2 -d)(rf- 4)4 (45) 



Tr[HG{r)f^d^ (46) 
r* 

These examples demonstrate that the asymptotic behavior of the correlations of 
critical points depend on the type of density, but are also sensitive to the particular 
value of 7 and the dimensionality d. For instance, for power law correlations with 
7 = 1, the slower decaying term (|^^ vanishes at d = 3, while for the logarithmic 
correlation vanishes at d ^ 2. 
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Appendix A. Details of the derivation in Sec. 2 



In this section we provide the details of the derivation of the results brought in Sec. 
[2] To begin with, we present the form of the matrices introduced in Sec. 12. ![ Eqs. 
(jlOP and As a first step, redefine the gaussian vector from Eq. ([5]): 

* = (V0(ri),V0(r2),t/,i,V2)- (A.l) 

This reordering of elements is more convenient in terms of the algebra to be carried 
out shortly. Then we have: 

A M 

a(0) air) \ f ^ ^(0 ^ \ 
air) a(0) ) \ -m(r) 
-TO(r)* \ / 6(0) b{r) 

V V '"W* / V ^(0) J J 



B = 



/ 



(A.2) 



(A.3) 



with the definitions 

a{r)ij 



b{r) 



ij.kl 



dridrj 
d'G{\r\) 



r=(r,0---,0) 



m(r) 



i,jk 



dridrjdrkdri 
d^G{\v\) 



= (r,0---,0) 



dridr.drk 



(A.4) 
(A.5) 
(A.6) 



= (r,0---,0) 

Here, r = r2 — ri, and b and m have double indices to denote elements of xp^. To make 
things clearer, we show some of these matrices explicitly for the three-dimensional 
case: 



a(0) = h G(2) (0) 



where Id is the three dimensional identity matrix, and 



6(0) 




r=(r,0---,0) 



= (r,0---,0) 

It can also be verified that a{r) and 6(r) have non-zero elements exactly for the same 
indices as a(0) and 6(0). We plug these definitions into Eq. and perform the 

integration over the field gradients and their Fourier conjugates. This leads to Eq. 
(dH) with 

B = I3-D (A.7) 

■ di(r) d2{r) " 
d2{r) rfi(r) 



D = APA-~H.I = 



(A.8) 
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where 

d,ir) = [a(0)2 - a(r)2] (-l)^-^a(s) (A.9) 

i = l 
r i = 2 



Long range asymptotics 

In order to obtain the long range asymptotics to second order in 77 (see Eq. ([T9|) ). we 
expand the prefactor and gaussian exponent in Eq. ()14p : 



1 



tvBtv 

2 



exp 

detA-1/2 



{(deta(O))^ 
1 



77 • -tzABtz 

' 2 



1 



rAB7 



,1 



-tvDtv 



■ dct 



/- 



(a(O)-ia(r))']} 



'1/2 



1 + itr (a(O)-Vr))' 



AS = 
D ~ 



deta(O) 
where here we have: 

b{r) 
b{r) 

di{r) 
di(r) 

since at the long range limit d2 ^ di, as follows from Eq. (|A.9p . Then we obtain Ci 
from Eq. ((22| and C2 from Eq. (|25p . which we can now write in full: 



C2(r) 



Co Tr [a{0)-'a{r)] ' - Cq Tr(Boi) + (p)"' Tr(FL) 



Co Ti{BoKBoK) - 2 Tr(FA'BoA') + Cg^i Tt{FKFK) (A.IO) 



with L = B^^DB^^ and A', F as in Eqs. (gS]) and ([Ml)- We note however that if G{r) 
diverges at the limit r — >■ 00, then the detA^^^^ prefactor must be expanded to order 
rf so as to get the proper subleading behaviour of C(r). 
Finally, let us examine F. We can write it as: 

(p(ri)t/'iV'i>i(p(i-2))2 (p(ri)V'i)i(p(r2)i/'2)2 
(p(ri)t/ji)i(p(r2)-f/'2>2 (p(ri))i(p(r2)tA2V'2)2 



F cx 



(A.ll) 



where 



I exp 



(A.12) 



V^det [27r6(0)] 

Now, recall that p cx dctiJ. From symmetry properties of the determinant, one can 
verify that the only non-zero in {ptp) (the position no longer matters) is 

1 



and the only non-zeros in (pxpxp) are (i ^ j): 
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To obtain the short range asymptotics, we assume G(r) is an even function of r and 
expand it to sixth order: 



1 2 1 4 1 6 

792r + —.gir - —g^r 



2r 4!^ 6!' 

where —(72,34, and —Qq are respectively the second, fourth and sixth derivatives of 
G{Q). Next, we need to find the prefactor in Eq. (|14p and transform B to an 
appropriate form. To find the prefactor we expand det A at the hmit r ^ 1. This 
yields: 

n2N , (52.94?-^) , 



AetA= \[{a{Q)l~a{rt 



After performing the rotation 
dimensions: 



an explicit calculation gives for three 
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where 60 = &i + A5. It is now useful to diagonalize 60 ■ Let us rewrite the vector 
as defined in Eq. (|29|). in the form 

■0+ = ■011 +0^, 

where 

0|| = (*ii,$i2,...,$ld,0,...), 
0^ ^ (0,...,$22,$23,---,$d<i)- 

Here -i/jy holds the first d elements of tjjj^, or in other words all those elements involving 
derivatives along the line connecting ri and Y2 ■ ip± holds all the other elements namely 
those involving derivatives in direction transverse to the vector r2 — ri. We then write 
the cross correlation matrix of xp, in the form 



bo = 



Co 



Cl 
C2 



where cq = (t/'y (E)ipi\),ci = (0|[ ® ipj 
be block-diagonalized: 



and C2 = (E) This matrix can now 



Co 



C1C2 ^c* 







C2 



C3 
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with 

/ 

CiC^"^ I 

Although in principal this requires transforming to new coordinates, in 

practice, cq = 0(7''*),ci = 0(r^),C2 = 0(1) so that to fourth order in the r expansion 
of C(r), u can be taken as the identity matrix. 

Substituting these expressions for the vectors and matrices into Eq. pi]) , 
performing the inverse Fourier transform, and rescaling the coordinates as: 

xp_ rxp_, Ab r^'^Ab 
converts B into a block-diagonal form, and thus 

This form is convenient for expansion in r, yet in order to evaluate the integral Eq. 
(|14p one should also expand p as a power series in r. Here care must be taken, as p is 
not a smooth function. Let us define a shorthand notation: 

MV^)=dcti/[</)(r,)]. 

so that 

where Q is some function which depends on the type of the density we are dealing 
with (unsigned or minima). Then we can expand the Hessian determinants as 

,,,,,,,, , ,'tp± + r^i'w + rtp^ \ fip± + r^ip^^ - rip_ \ 
h (V'l) h {ip.^) h h 




V2 J \ V2 J 

^) + ■ [hV ®Vh-Vh® Vh] (fy \-ip_+0 {r^^ 



Now, recall that the Hessian matrix one can construct from the elements of is 















$22 


$23 





$23 


$33 



(A.13) 

This implies that h (^"^^ = since both the first row and first column equal zero. In 
addition, of all the elements of the vector V/i only the first element survives: 

dot H = det Hi_ = det ( J^^ ) (A.14) 

9*11 V *23 $33 y ^ ' 

This formula holds for any dimension. It allows us to write for the unsigned 
correlation: 



(p(O)p(r)) 



2 

-y V- (V/i (g) Vh) ip_+0 (r*] 



-Afen^(dctFi)') + 0(r'*) (A.15) 
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thus recovering our main result, Eq. with: 



as 



3— 



-•Afeii((detff_L 



(A.16) 



Note that ^(det_ffj^) y is the expectation of a polynomial in gaussian variables and 
can be evaluated with Wick's theorem. 

For the minima-minima correlation we find: 

^2 



[p{0)pir) 



X e 



-■ijj_{Vh®Vh)'4>_ + O I 



d-2 



]Je (dcti7_L 



i=Q 




O 



To obtain this result we have used the fact that ijj_{'Vh (E) 'Vh)ip_ is negative 
semidefinite. Also, note the fact that all the minors of H at both points can be 
approximated by H± so that the two sets of 8-functions (from the two densities) 
merge. 

Obtaining the next order is considerably harder. However, we can deduce its 
power law dependence from the following considerations. Let us concentrate on the 
minima-minima correlation, so that: 

{[p (0) p (r)]) ^ r' { [-x'hl + ^r^] 6 [-x'hl + pr'] ) 

Here x is the first element of h± = dct H±, and p is some function of i/jj^, 
and xp_. Now, for the O function to be nonzero we require 



^<r^«l. 



(A.17) 



We can essentially ignore p in this equation since as long as p is not of order ~ 1 
this condition doesn't hold, so we have x^Va. -^i r^ <^ 1. Let us change coordinates by 
replacing x with x = x?7_l, keeping all other variables. This gives: 



([p(0)p(r) 



dx{{ 



- pr^] e [r 



r' I ^ + ^ 



(A.18) 



and therefore C (r) cx r^~'^ which is the result stated in Eq. ([S]) for fc = 3. For 
the unsigned correlation, there is no condition such as (jA.17p . so the integration will 
simply give a correction of order r^ above the first order. This can be seen by noting 
the identity |a;| = 2x Q {x) — x. 



Application of short order analysis in two and three dimensions 

Let us apply the short-distance result expressed in Eq. (|A.15p to the two and three 
dimensional problem. To this end we apply Wick's theorem to Eq. (|A.16p and find 
that, in two dimensions, the unsigned correlation behaves as 

C (r) ^ as + O {r^) (A.19) 
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with 

«3 - [kl - kl) , (A.20) 

where /cq = 54/521 s^nd = 55/521 while for the minima-minima correlation we have 
C (r) oc r^. 

In the three dimensional case the short range asymptotics of the unsigned 
correlation function takes the form 

C(r) = 03- + 0(r) (A.21) 
r 

with 

29 



«3 = ^v/fco(fc?-fco'), (A.22) 

while for the minima-minima correlation we find C(r) oc r^. 

To verify our analysis, we compare our results to numerical calculations. The 
correlation functions chosen are the same as those appearing in sections 13.11 and 13.21 
Figure lAll shows the excellent agreement we obtained with numerical calculations. 
There is one problematic graph, namely lAlbl The reason for this is simply the 
technical difficulty in numerically extracting the next-leading order behavior for this 
problem. 

Appendix B. Details of numerical analysis 

In this section we outline the methods we used to evaluate correlations numerically. 
From Eq. (fT4| in Sec. 12. H it follows that one can represent the correlation in the form: 



C(r) = ,,^r„ ^ , ^ / #p(0)p(r) exp f-it/.^- V ) (B.l) 



(27r)''VdetA (27r)2("-d) Vdct S 



2 



Here the elements of both A and B are functions of r. Evaluating the det A term 
numerically poses no difficulty. In order to evaluate the integral, we perform a 
Cholesky decomposition [5S] numerically on B: 

B = yV (B.2) 
tl^^V-^ip (B.3) 

This transforms the integral to an integral over an appropriate set of standard normal 
variables: ^ N{0, 1). The integral can then be estimated by standard Monte-Carlo 
integration. Our integrations were performed with Sobol numbers of the appropriate 
dimension generated by Mathematica 7. 

The evaluations of the two-dimensional problem required 2^^ lattice points for 
every data point appearing in figure [1] and 2^* lattice points for every data point in 
figure lAll For the three dimensional problem, we used 2^^ lattice points for data 
points in figure[51 and between 0.875 and 1.125 x 2^^ points (depending on the specific 
dataset) for data points in figure ED 
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(c) Minima-minima correlation in 3D. (d) Unsigned correlation in 3D. 
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